Transmission resonances for a Dirac particle in a one-dimensional Hulthen potential 
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We have solved exactly the two-component Dirac equation in the presence of a spatially one- 
dimensional Hulthen potential, and presented the Dirac spinors of scattering states in terms of 
hypergeometric functions. We have calculated the reflection and transmission coefficients by the 
matching conditions on the wavefunctions, and investigated the condition for the existence of trans- 
mission resonances. Furthermore, we have demonstrated how the transmission-resonance condition 
depends on the shape of the potential. 
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The study of low- momentum scattering in the Schroinger equation in one-dimensional even potentials shows that, 
as momentum goes to zero, the reflection coefficient goes to unity unless the potential V(x) supports a zero-energy 
resonance In this case the transmission coefficient goes to unity, becoming a transmission resonance 0- Recently, 
this result has been generalized to the Dirac equation[3|, showing that transmission resonances at k = in the Dirac 
equation take place for a potential barrier V = V(x) when the corresponding potential well V — —V(x) supports 
a supercritical state. The conclusion is demonstrated in both special examples as square potential and Gaussian 
potential, where the phenomenon of transmission resonance is exhibited clearly in Dirac spinors in the appropriate 
shapes and strengths of the potentials. Except for the both special examples, the transmission resonance is also 
investigated in the realistic physical system. In Ref.|4|], a key potential in nuclear physics is introduced, and the 
O r scattering and bound states are obtained by solving the Dirac equation in the presence of Woods-Saxon potential, 
which has been extensively discussed in the literature [f| H, @, B H|- The transmission resonance is shown appearing 
at the spinor wave solutions with a functional dependence on the shape and strength of the potential. The presence 
C$ . of transmission resonance in relativistic scalar wave equations in the potential is also investigated by solving the 
one-dimensional Klein-Gordon equation. The phenomenon of resonance appearing in Dirac equation is reproduced at 
'— 1 , the one-dimensional scalar wave solutions with a functional dependence on the shape and strength of the potential 
similar to those obtained for the Dirac equationfioj. Recently, Villalba and etal. have discussed the scattering of 
a relativistic particle by a cusp potential [ll|. They have solved the two-component Dirac equation in the presence 
s^q of a spatially one-dimensional symmetric cusp potential, and derived the conditions for transmission resonances as 
f-^. \ well as for super- criticality. Similarly, they have also solved the Klein-Gordon equation in the presence of a spatially 
I/"") ■ one-dimensional cusp potentialjl^], and obtained the scattering solutions in terms of Whittaker functions together 
' with the condition for the existence of transmission resonances. 

Due to the transmission resonance appearing in the realistic physical system for not only Dirac particle but also 
i— I ' Klein Gordon particle as illustrated in the Woods-Saxon potential as well as the cusp potential, it is indispensable to 
C***~ . check the existence of the phenomenon in some other fields. Considering that the Hulthen potential [l3| is an important 
realistic model, it has been widely used in a number of areas such as nuclear and particle physics, atomic physics, 
condensed matter and chemical physics[I3, [H, HE 03 ■ Hence, to discuss the scattering problem for a relativistic 
particle moving in the potential is significant, which may provide more knowledge on the transmission resonance. 
Recently, there have been a great deal of works to be put to the Hulthen potential in order to obtain the bound and 
scattering solutions in the case of relativity and non relativity [l8j. However, the transmission resonance is not still 
checked for particle moving in the potential in the relativistic case. In this paper, we will derive the scattering solution 
of the Dirac equation in the presence of the general Hulthen potential, and show the phenomenon of transmission 
resonance as well as its relation to the parameters of the potential. 

According to the definition in Ref. [HI, [l!| , the general Hulthen potential is chosen as 

v{x) = Q{-x)^— + e( x )^—, (i) 

where all the parameters Vb, a, and q are real and positive together with q < I being required to remove off the 
divergence of Hulthen potential, i.e., < q < I here. If q = 0, the Hulthen potential will degenerate to the cusp 
potential represented in Ref.pjJ, where the transmission resonances have been discussed in details. 0(x) is the 
Heaviside step function. 
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In order to investigate how the transmission resonance happens in the potential for a Dirac particle, following a 
similar procedure to that used by the Dirac equation takes the form (h = c = 1) 



"f— - 

\dx» 



iqA,_ 



ml 



* = 0, 



(2) 



where the four-vector potential can be written in a covariant way as qA^ = iV{x)§J t with 7 4 = (3 and x 4 = it 
there. When limited to the 1+1 dimensions, the Dirac equation (1) in the presence of the spatially dependent electric 
field becomes 



l 1 -^-i!3—+f3V{x)+ml 



# = 0, 



(3) 



In Eq.(3), V(x) is independent on time, hence the time dependence of the spinor "J can be separated with \P = e iEt ip 
as follows 



7 1 -^- - (E - V(x)) /3 + ml 
ax 



ip = 0. 



(4) 



Taking into account that we are working in the 1+1 dimensions, it is possible to choose the following representation 
of the Dirac matrices, i.e., taking the gamma matrices 7 1 and (3 to be the Pauli matrices a x and a z respectively, then 
the Dirac equation truns into 



o x - (E - V(x)) a z + ml 

ax 



ip{x) = 0. 



The four-spinor, ip(x), is decomposed into two spinors, u\ and U2, so that 

ip(x) ■■ 



ui (x) 
u 2 (x) 

Thus the problem is to solve the coupled differential equations: 

u[(x) = - (m + E -V(x))u 2 (x) , 
u 2 (x) = — (to — E + V{x)) u\ (x) . 

By introducing the following combinations 

4> (x) = ui (x) + iu 2 (x) , x ( x ) = u i ( x ) - i u 2 (x) ■ 

Substituting these into (7) and (8) and re-arranging gives: 

(f)' (x) — —imx (x) + i(E — ^(a:)) (j> (x) , 
X' (x) = imcf) (x) — i (E — V (x)) x (x) ■ 

The two components, (j){x) and x ( x ), satisfy: 

0" + \{E - V) 2 - m 2 + iV' (j) = 0, 
x" + \{E - V) 2 - to 2 - iV'} x = 0. 



(5) 



(6) 



(7) 



(9) 



(10) 
(11) 



(12) 
(13) 



In order to obtain the scattering solutions for x < with E 2 > 1, substituting the potential presented in (1) into 
(12) gives 



d 2 4>{x) 
dx 2 



E- 



V 



a V Q ae~ a * 
— to + 1- 



(e~ ax - q) 



(x) = 0. 



On making the substitution y = qe ax , Eq.(14) becomes 



2„,2 ' 



a y 



$ + a 2 — 
dy 2 dy 



E 



Vo y 
q 1 - v 



iVpa V 

1 (i-yf 



(x) = 0. 



(14) 



(15) 



In order to derive the solution of Eq.(15), we put <f> = y^ (1 — y) /, then Eq.(15) reduces to the hypergeometric 
equation 



V (1 - V) f" + [1 + 2 M - (2/i + 2A + 1) y] !' - (^A + ^§^j f = 0, 

where the primes denote derivatives with respect to y and the following abbreviations have been used 

ik ip iVo 

H = — , v = — , A = , 

a a aq 

p 2 = (E + V /q) 2 -m 2 , k 2 = E 2 -m 2 . 



(16) 



(17) 



Note that as we are considering scattering states, \E\ > m which ensures that k is real, and Vq is real and positive, p 
is real for q > 0. The general solution of Eq.(16) can be expressed in terms of hypergeometric function as 



f(y) = A F (n - v + A, /x + v + A, 1 + 2fi; y) + B y'^F (-fi - v + A, -fj, + v + A, 1 - 2fi; y) . 



(18) 



So 



b L {y) = A y» (1 - yf F (fi - v + A, ^ + v + A, 1 + 2 M ; y) + B y~» (1 - y) X F (-p - v + A, - M + v + A, 1 - 2^; y) . 



As x — ► — oo, there is y — > 0. So, the asymptotic behavior of 0l(j/) can be written as 

lim 4> L (x) = Aq lk / a e lkx + Bq- ik ' a e~ ikx . 



x > — OO 



(19) 
(20) 



From equation (10) the other component, x( x ) 1S 

1 



X (x) = ^-[i (E-V(x))cj>(x)-cf>'(x)} 



(21) 



Substituting equation (20) into the above gives us 

.:e-i 

hm xl (x) = A 



x > — oo 



ik/a £ ikx _|_ q 



E + k 



—ik/a —ikx 



(22) 



The choice of combinations of the wave function components (9) can be re-written : 



(23) 



Upon substitution of equations (20) and (22) into the above it can be seen that the wave function, ip(x), comprises 
of an incident and reflected wave far to the left of the barrier which is the desired form to establish reflection and 
transmission amplitudes. 

Next, we consider the solution of Eq.(3) for x > 0. With the potential represented in Eq.(l), the differential 
equation to solve becomes 



d 2 4>{x) 
dx 2 



+ 



E 



Vn 



2 • V L 

m — i- 



(e ax - q) 



(j)(x) = 0. 



(24) 



The analysis of the solution can be simplified making the substitution z = qe ax . Eq.(24) can be written as 



2 2 d 2 <f> 2 d<t> 
a^z^—^ + a z — 
dz z dz 



E 



Vo z 
q 1 — z 



m 



iVoa z 
~(1-^) 2 



<j)(x) = 0. 



(25) 



Put = z^ (1 — z) g, Eq.(25) reduces to the hypergeometric equation 



:(l-z) g" + [1 + 2/i - {2pi - 2A + 1) z] g' - 



a 2 q J 



.9 = 0, 



(26) 
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where the primes denote derivatives with respect to z. The general solution of Eq.(26) is 

g(z) = C (jjl - v - A, fi + v - A, 1 + 2^; z) + Dz~ 2fl F (-/z - v - A, -fi + v - A, 1 - 2/x; z) . (27) 

So, 

(j) B (z) = Cz» (1 - z)~ X F(ji-u-X,fi + u-X,l + 2/i; z) + Dz~^ (1 - z)~ A F (-fi - v - A, -fi + v - A, 1 - 2/z; z) . 

(28) 

Keeping only the solution for the transmitted wave, C = in Eq.(28). As x — > +oo(z — > 0), there is 

fo(a;) — ► Dq~ ik ' a e ikx . (29) 

while the other component 

Xfl(as) — > D f— q- ik / a e ikx . (30) 
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The electrical current density for the one-dimensional Dirac equation is defined as 

j =$(x)~f x il)(x) = -ipix)^ a 2 ip(x) =i(u*u 2 - u^ui) = i (|0(a;)| 2 - |x(^)| 2 ) (31) 

The current as x — ► — oo can be decomposed as = j m — j Ic n where j- ln is the incident current and j re fl is the 
reflected one. Analogously we have that, on the right side, as x — ► oo the current is ju = jtrans, where jt ran s is the 
transmitted current. Using the reflected j re fl and transmitted Jtrans currents, we have that the reflection coefficient 
R, and the transmission coefficient T can be expressed in terms of the coefficients A, B, and D as 

jin \A\ 2 \E-kJ v 1 

T = ^ = (33) 

jin \A\ 2 

Obviously, R and T are not independent; they are related via the unitarity condition 

R + T=l. (34) 

In order to obtain R and T we proceed to equate at x = the right 4> R and left <j> L wave functions and their first 
derivatives. From the matching condition we can derive the relations among the coefficients A, B, and D. Then the 
reflection coefficient R and transmission coefficient T are obtained. 

The calculated transmission coefficient T varying with the energy E is displayed in Figs. 1-4 at the different values 
of the parameters in the Hulthen potential. From Figs. 1-4, one can see that the transmission resonance appears in 
all the Hulthen potential considered here. But the intensity and width of resonance as well as the condition for the 
existence of resonance are different, and they depend on the shape of the potential. Compared Fig.l with Fig. 2, 
it can be seen that the width of resonance decreases as the decreasing of diffuseness a, which is similar to that of 
Woods-Saxon potential as shown in Figs. 3 and 5 in Ref . [Io| . The same dependence can also be observed from Figs. 3 
and 4. Compared Fig.l with Fig. 3, one can find that the condition for the existence of transmission resonance does 
also relate to the parameter q. As q decreases, the height of potential barrier increases, the widths of the transmission 
resonance increases. The conclusion can also be obtained by comparing Fig. 2 with Fig. 4. In order to obtain more 
knowledge on the dependence of transmission resonance on the shapes of the potential, the transmission coefficient 
T varying with the strength of potential Vq is plotted in Figs. 5 and 6. Beside of the phenomenon of transmission 
resonance, similar to the Fig.l and 2, the width of resonance decreasing as the decreasing of diffuseness a is disclosed. 
All these show the transmission resonances in Hulthen potential for Dirac particle possess the same rich structure 
with that we observe in Woods-Saxon potential. 
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FIG. 1: The transmission coefficient for the relativistic Hulthen potential barrier. The plot illustrate T for varying energy E 
with Vo = 4, a = 1, and q = 0.9. 
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FIG. 2: Similar to Fig.l, but with V = 4, o = 0.5, and q = 0.9. 
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FIG. 5: The transmission coefficient for the relativistic Hulthen potential barrier. The plot illustrate T for varying barrier 
height V with E = 2, a = 1, and q = 0.9. 



